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1. Introduction 

The low dimensional fermionic models as toy models provide an excellent playground to try 
and develop new ideas and methods in quantum field theory jy]. These models have only 
indirect physical meaning but they are much simpler than their 4-dimensional counterparts, 
and they usually show important characteristics of the original ones. For example the 2- 
dimensional quantum electrodynamics (QED2) proved to be a good toy model for treating 
the soft mechanism of the quark confinement |2|, j3|. The confining properties, the large- N c 
expansion Q or the baryon structure [||, ||] of QCD can also be studied in the 2-dimensional 
version of the model and then one can get even analytical results for the non-perturbative 
domain. One usually takes the bosonized version of these models which are local self- 
interacting scalar theories, and can be investigated in an easier way J?]. 

The phase structure of the QED2 with many flavors was mapped out from its bosonized 
version and it was shown that it exhibits only a single phase || ||] as opposed to the 
single-flavor QED2 (which is often referred to as the massive Schwinger model) M, [3|], 
which possesses a symmetric strong coupling (e 3> m e ) phase and the weak coupling 
(e <C m e ) phase with spontaneously broken reflection symmetry separated by the critical 
value (m e /e) c ~ 0.31 as was shown by density matrix renormalization group (RG) technique 
1 1C ] or by continuous RG method [|ll| ^ . 

The situation is a bit different in the case of the 2-dimensional quantum chromody- 
namics (QCD2). The scalar model equivalent of the single flavor QCD2 can be easily 
obtained, but the bosonization failed to treat the multi-flavor model. The solution of this 
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problem was the introduction of the non-abelian bosonization [12]. However a single- flavor 
but multi-color investigation is possible with the usual abelian bosonization technique. 

It is argued [13] that the multi-color QCD2 possesses two phases, a weak coupling or 
quark phase, and a strong coupling or Bose phase. Elsewhere it is also argued Q] that 
the model has a single phase. This open question should be investigated in the low energy 
limit. 

The bosonized version of the multi-flavor QED2 contains sine-Gordon (SG) type pe- 
riodic self- inter action terms. The scalar fields are coupled by a mass matrix giving a 
multi-component or layered sine-Gordon (LSG) model which is used to describe the vor- 
tex dynamics of magnetically coupled layered superconductors [15], where the number of 
flavors in QED2 equals the number of layers of the condensed matter system II. The Bose 
form of the multi-color QCD2 also contains SG type interactions, a mass matrix and a 
mixed term. The latter can be associated to the non-periodic term in the potential there- 
fore it can be Taylor expanded giving further contributions to the mass matrix. The higher 
order terms are negligible since they do not modify the phase structure of the model. Then 
one can conclude that the difference between the bosonized versions of QED2 and QCD2, 
respectively stems from the different mass matrices. For the 2-flavors (Nf = 2) QED2 and 
the 2-colors (N c = 2) QCD2 the mass matrices coincide, implying that these models are 
equivalent. It is quite surprising since the fermionic models are different. The difference 
between the bosonized models appears when Nf > 2 and N c > 2. We note that the low 
energy multi-flavor QCD2 with unequal masses can also be bosonized giving a so called 
generalized SG model || [R| but its investigation is out of the scope of this paper. 

The phase transition of these models was obtained from the microscopic theory so far, 
which is formulated in the high energy /ultraviolet (UV) region. The low energy /infrared 
(IR) physics can be obtained by integrating out the quantum fluctuations one by one. Then 
one can get a low energy theory describing the quantum system at energy scales where 
the measurements are usually performed. The quantum fluctuations can be eliminated 
systematically by using the renormalization group (RG) method. The original fermionic 
models contain strong couplings in the high energy UV regime which disables one to perform 
a perturbative renormalization. The evolution is usually started from a perturbative region 
where the theory is almost interaction free. We also note that the RG equations should 
preserve the gauge symmetry [17]. However the bosonized version of the toy models, which 
are simple scalar models can be easily treated by the functional RG method. The low 
energy IR physics of the 2-dimensional one-component scalar field theories which contain 



periodic self-interaction term are well understood 11, 18, 19]. 

Our goal in this article is to compare the phase structure of the multi-flavor QED2 
and the multi-color QCD2 by using functional RG method. The bosonization is applicable 
for the original fermionic models at a certain parameter choice /3 2 = 4-7T, which appears in 
the argument of the sine function. Therefore the value of the parameter /3 has to be kept 
fixed during the RG evolution. It implies that the investigation which is confined to the 
local potential approximation (LPA) where /3 does not evolve may give reliable evolutions. 
Although the Wilsonian renormalization procedure is very powerful in LPA [ 20 ] , we choose 
the effective average action RG method f2l| , 22, 2^, 19] to obtain the evolution of the models 
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due to its more flexible usage. The RG evolution provides us the scale dependence of the 
couplings in the scalar model from which the scaling of the original fermionic couplings can 
be obtained according to the bosonization rules. 

The paper is organized as follows. In Sect. ^| we introduce the bosonized versions of 
the investigated fermionic models and relate them to layered sine-Gordon (LSG) models. 
The evolution of the couplings is determined in Sect.|3| Finally, in Sect. |] the conclusions 
are drawn up. 



2. Bosonized models 



2.1 Multi-flavor QED 2 

The Lagrangian of the multi-flavor QED2 is 



N f N f 



(2.1) 



i=i 



i=l 



with Nf Dirac fields and identical fermionic charge e and mass m e , furthermore Fq\ = 
OqAi — OiAq. One can transform the fermionic field variables ipi,ipi into bosonic ones (pj 
by the bosonization rules [24, 25, 0] 
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sin (2^^), 



7T 



(2.2) 



where N me means normal ordering with respect to the fermion mass m e and c = exp(7)/27r, 
with the Euler constant 7 = 0.5774. The Hamiltonian of the system in Coulomb gauge is 
given by 

A ^ 



Nf 2 
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1 J x ^ J x 
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with f = Jq dx° J^ L dx l and 
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The resulting bosonized form of the Hamiltonian is 



[x . 



(2.4) 



U = J\frr 
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where Ilj(x) denotes the momentum variable canonically conjugated to 4>i(x). In or- 
der to complete the bosonization, one has to use normal-ordering with respect to the 
scalar mass /i 2 = e 2 /-7r which modifies the coupling of the periodic term, M me cos ( 
(^/m e ) fc2//47r A/' At cos(6</>). Therefore, the Nf = 1 flavor bosonized QED2 reads as 



/'I 1 1 

-n 2 (x) + -{dict){x)) 2 + -n 2 (f>(x) 2 - cm e n cos (2^4>{x)) 



(2.6) 



which can be generalized for Nf > 1 flavor using its rotated form where the mass matrix 
is diagonal. Let us note that the scalar mass term in (^]5|) can be rewritten in terms of the 
mass matrix -Mq ED defined via 

2 f N f \ 2 

-$A^ ED $ = -f- f^a n ^j , (2.7) 

where <3? = (</>i,</>2) ■ ■■,4'Nf) and the couplings a n are free real parameters of the model. In 
order to reproduce the mass term of ( |2.5| ) one has to restrict the choice to be a n = 1 for all 
n. However, based on symmetry considerations any choice with a 2 = 1 for all n = 1, . . . , Nf 
should reproduce exactly the same phase structure since the number of zero and non-zero 
eigenvalues of the mass-matrix remains unchanged which was found to be decisive with 
respect to the phase structure of the iV/-component model [[|. It is not a surprise since 
in path-integral quantization a change of sign of any of the field components represents an 
allowed transformation of the integration variable, that in turn induces a single change of 
sign in the mass term and leaves the other terms of the action invariant. As a consequence, 
a n = (— l) n+1 is also a suitable choice which will be used in the mass term of QED2 in this 
paper and leads to the mass matrix 

(•MqedW = (-l) a+b G, a, b = 1, 2, . . . , N f (2.8) 

with G = e 2 /vr. The QED-type mass matrix ( |2.8[ ) exhibits a single non-vanishing mass- 
eigenvalue, Mjy- = NfG and Nf — 1 vanishing eigenvalues. 

2.2 Multi-color QCD 2 

The Hamiltonian of the QCD2 with a single flavor Nf = 1 is 

N c Nc N c 

■H = g 2 Y^E b a 2 +Yl ^71(^1 - A b a )i> b + m g ^2 ^Va (2-9) 

a, 6=1 a, 6=1 a=l 

in the gauge 

A = 0, A a b = for a = b, El = for a ^ b. (2.10) 
Using the Gauss law the bosonized Hamiltonian with one flavor becomes 



I (13 + (<W„) 2 ) - ^A/-„co S (2v5A,) 
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a,6 a, 6 



- 4 - 



where the scale [i should satisfy ^ = c'g, with d a constant, in order to take the interaction 
energy proportional to g 2 fi~3[| . 

We determined the field configuration for the ground state of the model numerically. 
The static field configuration minimizing the energy was searched for by means of the con- 
jugate gradient method as in ||. The results showed that the constant field configuration 
with all components <j) a equal to the same constant minimizes the energy, for the cases 
N c = 2, 3. The same trivial ground state appears when we take the polynomial piece of 
the potential alone. Therefore, it seems to be justified to Taylor-expand the last potential 
term in the Hamiltonian at <j) a — </>;> = and keep only the quadratic term of the expansion, 
so far one is interested in the behavior of the system in or close to the ground state. Then 
one finds 



H 



1 



nl + (Ma 



CC 



m g g 



1 



Mfj, cos(2^4> a ) + ^ g 2 c g (4>a - 4>bY 

b 



(2.12) 



with c g = 1/AitN c + (4cc') /3. Let us note, that similarly to the bosonized QED2, the 
scalar mass term in ( [2.12|) can be rewritten in terms of the mass matrix A^q CD defined via 



QCD 



a,b 



(2.13) 



where 



!>i, (f>2, •••) 4>N C ) and the summation runs from a, b = 1 to N c . Then one gets 

(M 2 QCD ) a M = (N-l)J5 a , b -J a,b = 1,2,..., N c (2.14) 



with J = 2g c g . The QCD-type mass matrix (2.14) has a single zero eigenvalue and N c — 1 



identical, non-vanishing eigenvalues, Mjy- = N C J. 
2.3 Relation to layered sine-Gordon models 

Both models, the bosonized multi-flavor QED2 and the multi-color QCD2 can be considered 
as the specific forms of a generalized LSG model [Hj which consists of two-dimensional 
periodic scalar fields coupled by an appropriate mass matrix whose bare Euclidean action 
is written as 



r\i 1 N 

^ x . n=l 



(2.15) 



with the O(N) multiplet $ = . . . , 0j\r)- F° r t ne specific choice (3 2 = An with the 
mass matrices ( |2.8[ ) and (|2.14| ) one one recovers the bosonized version of the multi-flavor 
QED2 for N = Nf and that of the multi-color QCD2 for N = N c , respectively. The 
amplitude y of the periodic piece of the potential is identical for all component fields, and 
it is proportional to the fermion mass (y ~ m), the exact relation can be determined by 
using normal-ordering w.r.t. the boson mass. We note that for N = 2 the mass matrices 
(|2.8| ) and ( [2.14 ) coincide, consequently, the Bose forms of the two-flavor QED2 and the 
two-color QCD2 are the same. 

For later use it is worthwhile mentioning that after an appropriate 0(N) rotation 
diagonalizing the mass matrix, the LSG model with the QED-type mass matrix ( |2.8| ) 
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exhibits a single massive field and N — 1 massless ones. On the contrary, the LSG model 
with the Q CD-type mass matrix ( |2.14| ) shows up a single massless field and N — 1 massive 
fields of identical masses after such an O(N) rotation. 



3. RG approach for multi-component models 



The systematic removal of the quantum fluctuations can be performed by the evolution 



equation for the effective action [21, 22, 23, 19] 



kd k T k 



2 



kd k R k 

r[ 2) + R k 



(3.1) 



(2) - 

where IjL refers to the second functional derivative matrix of the effective action and the 
trace TV stands for the integration over all momenta. The scale k starts from a large UV 
value A (which is typically set to oo during the calculations) and goes to zero. R k plays 
the role of the IR regulator function. For the suppression of the high-frequency modes one 
can choose the power-law type regulator 



Rk 



with the parameter 6 > 1. Here we choose 6=1, 



(3.2) 

which corresponds to the Callan- 
Symanzik RG scheme 23]. It is easy to see, that in d = 2 the chosen CS scheme is free of 
UV divergences and ultralocal, furthermore the evolution equations take a rather simple 
form. We note that in d = 2 the choices 6 = 1 and 6 = oo coincide [^, ^7]] in the LPA. 
The latter case corresponds to the sharp cutoff limit, which makes the functional form of 
the CS and the sharp cutoff (or Wegner-Houghton type) evolution equations similar. The 
effective action is expanded in powers of the derivative of the field, 



(3.3) 



with Vfc[$] the potential and Z&[$] the wave- function renormalization. The latter provides 
evolution to the parameter /?, even in the case of field-independent wave-function renor- 
malization 2Tfc[<l>] = z k , where 1 = 1/z [O]. The bosonization gives a constraint to the 
parameter (3 2 = Air. The running of Zfc[3>] influences the evolution significantly in the 
vicinity of the Coleman point at /3 2 = 8tt Jl9| , 28] but gives slight modifications around 
/3 2 = 47r 1 1 1 1 j therefore it is not supposed to affect the phase structure of the model. Thus 
we do not go beyond the LPA, and set = 1. One can derive the evolution equation 

1 



(2 + k d k ) Vfc($) 



An 



hi 



det <% + 



(3.4) 



with V* J 



dfcdffr.Vk for the dimensionless potential V k 
homogeneous field configurations. We make the general ansatz 



k 2 V k , where $ stands for 



1 N 
Vfe($) = -<5>M 2 {k)<5> + y{k) cos(/3 cf> n 



(3.5) 



n=l 
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for the dimensionless potential of the LSG type models under discussion, where y(k) = 
k~ 2 y(k). Inserting the ansatz ( |3.5| ) into Eq. (3.4), the right-hand side becomes periodic, 



while the left-hand side contains both periodic and non- periodic parts f2|, |(J ||, g| . The 
non-periodic part contains only mass terms, so that we obtain a RG flow equation for the 
dimensionless mass matrix 

(2 + kd k )M 2 (k) = 0, (3.6) 

giving the scaling 

j k = k- 2 J, and G k = k- 2 G, (3.7) 

which corresponds to the scaling according to the canonical dimensions, since in the LPA 
the anomalous dimension is zero. One can conclude that the dimensionful couplings J, G 
remain constant during the blocking. 

The RG flow avoids the singularities if we handle the evolution without any truncations 



1 21]. We note however that one should usually use some approximations or expansions in 
the RG equations in order to solve them. The bosonized QCD2 contains a single Fourier 
mode, however the RG equations generate the higher harmonics. Restricting ourselves 
to follow the evolution of the fundamental mode only may induce a strong truncation, 
implying that we should face the problem of poles of the evolution equation in Eq. (|3.4|). 
The proper choice of the IR regulator function may drive the evolution to reach the pole 
only in the k — > limit. This seems to be true even in the case of the sharp cutoff scheme, 



where one can draw up the quantum cenzorship conjecture [32, pSf] . However the poles 
- unless they appear as artifacts of the approximations and truncations, - have a great 
physical importance, because their existence can signal the spontaneously broken phase of 



the model [21, £4|]. In the symmetry broken phase of a single component scalar field the 



dynamical Maxwell cut makes the effective potential superuniversal, namely 

V k =om = (3.8) 

3.1 UV scaling 

The correct UV scaling can be obtained if we improve the results of the linearized approx- 
imation by taking into account corrections of the order O( J) for the QCD2 type and those 
of 0(G) for the QED 2 type case ]^9|, |, [jOj. This is achieved by linearizing the RG 



equation in the periodic piece of the blocked potential, 

(2 + kd k )V k = + 0(V 2 ), (3.9) 

where C and Fi(U k ) stand for the constant and linear pieces of the determinant 

det[c% + V*>] =C + iq(Vfe) + 0(V 2 ). (3.10) 
The UV scaling law for the QCD2 type LSG model for N colors is 



k\ N ^ ( k 2 + NJ 
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with the initial value y(A) at the UV cutoff k = A. From the extrapolation of the UV 
scaling law ( |3,11 ) towards the IR scales we can read off the critical value (3 2 (N) = 8irN. 



The coupling y is irrelevant for (3 2 > f3 2 (N) and relevant for /3 2 < /3 2 (N). The critical 
frequency and the corresponding critical temperature 

(iV ) _ 2tt 2. f o 12 N 

QCD ( AQ KTB jy (3-12) 

separating the two phases of the model coincide with the general expressions obtained 



previously for the rotated LSG model in Refs. [£9L 31, [3(J, 0L 

Similar consideration can be done for the QED 2 type LSG model ||, |l|, |, § and the 
solution of Eq. ( |3.9| ) for the couplings in case of A" flavors is given as 

/as/^A^ 3 " fk 2 + NG\ N ^ 
^) = "( A )(a) (iTSc) ■ < 313 > 

The critical frequency and the corresponding critical temperature which separates the two 
phases of the model can be read off directly, 

P^ N )- N _ V "> T QED " p 2 (N) ~ T KTB N ■ (3-14) 

For N = 1 layer the LSG model with magnetic type coupling reduces to the massive 2D- 
SG model where the periodicity is broken explicitly. Therefore, for sufficiently small bare 
coupling y(A), there exists only a single phase [^9|, |3(], [0| ^] in the N = 1 layer model, 
i.e., the coupling y(k) is relevant (increasing) in the IR limit {k — > 0) irrespectively of /3 2 . 
For — > oo the magnetically coupled LSG behaves like a massless 2D-SG model with the 
critical frequency /3 2 = 8tt ||. Both the dependencies in Eqs. ( 3.12|) and ( |3.14|) of the critical 



frequencies f3 c on the number of layers indicate that QED2 and QCD2 with any number of 
flavors and colors, respectively belong to the symmetry broken phase of the corresponding 
LSG model (f3 2 = Air < fil{N)), see also Fig.|. We shall show below that the IR scaling 
laws confirm this statement. The UV treatment suggests, that the phase structure of 
the LSG type models and the SG model is quite similar, since there is a Coleman fixed 
point with the critical parameter f3 c separating the symmetric (/3 > (3 C ) and the symmetry 
broken (/3 < f3 c ) phases. The Coleman point appears in the UV level, namely the couplings 
start to scale irrelevantly (relevantly) in the (broken) symmetric phases respectively. The 
inclusion of the higher harmonics modifies this picture by introducing further critical j3 



values [33], and may give further fixed points. However the IR behavior of the SG model 
results in a single Coleman point, suggesting the IR nature of the fixed point. 

Occurring a pole during the RG evolution may give another signal that the bosonized 
multi-flavor QED2 and multi-color QCD2 correspond to LSG models in the symmetry 
broken phase. The following search for such a pole relies on the extrapolation of the 
UV scaling laws again. Nevertheless, the appearance of a pole would inform one on the 
superuniversality of the effective potential settling it as ([^) with the consequence that 
the original multi-color QCD2 exhibits a single phase. 
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1 2 3 4 5 6 



N 

Figure 1: The critical frequency (3^(N) versus the layer- number N is shown for LSG models with 
QED2-type (dashed line) and QCD2-type (solid line) interlayer couplings, respectively. The critical 
frequencies lie outside of the shaded area, irrespectively of N. 



As stated earlier, QED2 for Nf = 2 and QCD2 for N c = 2 coincide. It was shown 
|8j that these models have a single phase. The situation changes for Nf = 3 and N c = 3, 
respectively. Then the RG equation in Eq. (|3.4[) for 3 layers has the form 



(2 + kd k )V k = ~ log 



;i + v k n )(i + vf )(i + v 33 ) + v k 12 v 23 v 31 + v k 13 v 21 v 32 



-V k 13 (l + V 22 )V 31 - V k 12 V 21 (l + Vf ) - (1 + V k n )V 23 V k 32 
For the 3-flavor QED2 the potential for the corresponding scalar model is 



(3.15) 



V k = -4>2 + hf + V [cosGS^i) + cos(/3<£ 2 ) + cos(/30 3 )] (3.16) 

with G > 0. Inserting this ansatz into the argument of the logarithm in the right hand 
side of the RG equation ( |3.15| ) one can see, that a pole for (f>\ = fa = 03 = may appear 
when 

(l-/3 2 y) 2 (l + 3G-/3 2 y) = 0. (3.17) 

Due to the factor in the first bracket in the left hand side, the relevant scaling of the 
dimensionless coupling y in Eq. ( p.!3| ) drives the flow to a pole independently of the UV 
initial parameters of the model. This simple treatment suggests that the model is in its 
symmetry broken phase, which implies that the 3-flavor QED2 has a single phase. Thorough 
calculations showed the same result |], j|] . 

The form of the potential for the 3-color QCD2 is 

Vk = - fa? + {fa - fa? + (fa - fa) 2 ] + y[cos(Pfa) + cos(f3fa) + cos((3fa)} (3.18) 
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with J > 0. The argument of the logarithm in Eq. ( 3.15 ) is now 

(l-/3 2 y)(l + 3J-/3 2 y) 2 = 0, 



(3.19) 



which also gives a pole due to the first bracket in the left hand side if we consider the 
relevant scaling of the coupling y according to Eq. ( |3.11 ). Similarly to the result obtained 
for the 3-flavor QED2 the model seems to be in the symmetry broken phase giving again 
a single phase for the original 3-color QCD2. Below we shall show that our expectation of 
a single phase for 3-color QCD2 is justified by the IR scaling laws. 



3.2 IR scaling 

A more reliable information on the phase structure of the multi-color QCD2 should be 
deduced from the IR scaling laws of the corresponding LSG type model. The proper 
treatment of the problem requires to find the solution of the rather complicated partial 
differential equation ( |3.4| ). Instead of treating that task in its full complexity, we shall 
invent the following strategy. First we perform an 0{N) rotation 1Z which diagonalizes the 
symmetric mass matrix by the rotated field variables a« = IZijfy. Due to the particular 
structure of the mass matrix -Mqcd' ^ exhibits a single zero eigenvalue and N — l identical 
eigenvalues NJ. Further on, we assume that the mass gap suppresses large amplitude 
quantum fluctuations of the massive field components and, therefore the potential can 
be Taylor-expanded in the massive field components at their vanishing value. Then the 
massive fields appear in the lowest order as free fields and decouple from the massless 
field component and can easily be integrated out. In this manner the problem becomes 
amenable for the numerical treatment. Such an approach has been successfully applied to 
LSG models in 

In particular, the rotation for N = 3 is performed with the matrix 



/ 



1 

1 

1 



1 

V3 



V3 



1 

75/ 



(3.20) 



and the effective action of the rotated model takes the form 



rot 



2^«i) 2 + \{d^2? + \{d^) 2 + -J{a\ + a|) + V rot 



(3.21) 



with 



t;, : OyooBl^jooB^jooB^ 



2^sin(^jcos^jsin^ 



+^cos ( ^) cos (/^) - ^cos (^) sin 
+y ~sin(^)cos(^)sin(^ 



(3.22) 
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Now we keep the lowest-order term of the Taylor-expansion of the potential in the massive 
field components at a<i = 03 = like in [31] for the LSG model. Then the effective action 
reduces to 

^red 



1 rot 



^{d^aif + 3ycos \ 



(3.23) 



The massive components decouple and describe free fields, and by integrating out the N — 1 
massive modes as described in Refs. || one obtains that the IR behavior of the model 
is completely determined by the remaining massless component field a% which describes a 
simple SG model. 

The advantage of the rotation reveals itself in the reduction of the task of the low-energy 
QCD2 to the determination of the IR behavior of the SG model ( 3.23[ ). The functional RG 



method showed ]35], |T^, 19] that this SG model has two phases depending on the value of 
its parameter /3 2 /3. As we have shown previously in Sect. |3.l| , the multi-color QCD2 gives 



values of the parameter (3 = Att < /3^(N) for arbitrary number N > 1 of colors, which 
means that the corresponding SG model is in the symmetry broken phase. Earlier calcula- 



tions based on Fourier expansion |35| , 36, 18, [26|] showed that in the symmetry broken phase 
the effective potential is superuniversal with the parabolic shape ([3.8]). It happened numer- 
ically that the Fourier-expansion drove the evolution towards the pole at a non-vanishing 
scale. Approaching it a parabolic prepotential appeared but the Fourier-expansion became 
unreliable at the same time, so that the further evolution was treated at tree level p7| 
which always gave parabolic effective potential ([3^). More precise calculations, avoiding 
any expansion of the potential J33| ] now suggest, however, that there is a non-trivial IR 
attractive fixed point for low values of j3 2 giving a superuniversal effective potential which 
deviates a little from Eq. (^^), and the latter form is reached only in the limit /3 2 — > 0. 

The flow of the coupling y is determined by a computer algebraic program [27], which 
solves the RG equation directly, without using any ansatz for the potential. It finds the 
fundamental mode y(k) by Fourier-analyzing the numerically determined potential at any 
scale k afterwards. The polynomial suppression scheme we use here needs higher numerical 
accuracy as compared to the exponential scheme pi]]. We set a high numerical working 
precision in order to handle the numerical ambiguities properly as was pointed in 33]. The 
IR flow of the coupling can be seen in Fig. ||. The figure clearly shows that the IR value of 
the dimensionless coupling y(0) is independent of its initial, microscopic value. This implies 
that the IR effective potential is also independent of the microscopic parameters, i.e. it 
is superuniversal. According to the inset of Fig. || the effective potential approaches the 
parabolic shape characterized by 1 + V"(3>) = in the limit k — > 0, but our high-precision 
calculation also shows that it does not reach the parabola, in accordance with the findings 
in ||. 

The dimensionful parameter J = g 2 c g is constant, therefore the fermionic coupling g 
remains unchanged during the evolution. It gives non-vanishing coupling in the IR limit 
of the 3-color QCD2. The IR behavior of the model shows that the dimensionful coupling 
y goes to zero for k — > driving the quark mass m g to zero. Therefore this scalar model 
is a free massive theory too as was the 3-flavor QED2 @, and the 3-color QCD2 is an 
interacting theory of massless two-dimensional quarks. 
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Figure 2: The IR scaling of the dimensionless coupling y for several initial values. The inset 
demonstrates the convexity of the effective action during the evolution as the scale k decreases, 
k/A = 10~ a . 



3.3 Large N case 

The techniques used in the previous subsections can be easily generalized to the case of 
arbitrary N. The naive expectation that multi-color QCD2 is in the symmetry broken 
phase stems from the appearance of poles if we consider the argument of the logarithm in 
the right hand side in Eq. ( |3.19 ) for N colors, 

(l-P 2 y)(l + NJ-p 2 y) N - 1 =0 (3.24) 

in the framework of the extrapolation of the UV scaling laws. Similarly to the 3-color 
case, the first factor can change sign due to the relevant UV scaling of the coupling y, 
signaling the appearance of the pole and the symmetry broken phase. However, a more 
reliable conclusion can be drawn again if one considers the IR scaling for iV-colors. In 
order to diagonalize the mass matrix, one performs the appropriate O(N) rotation. Taylor- 
expanding the potential in the new massive field variables a.{ , % = 2 . . . N at their vanishing 
values and keeping the quadratic terms only, these become massive free fields. Integrating 
them out one can reduce the effective action to that of the single massless field a±, 



■pred 
1 rot 



i(cy*i) 2 +Ny cos (p ° 



(3.25) 



This is a SG model too with decreasing parameter /?' = P/yN for increasing number N 
of colors. The calculation for smaller values of P' requires extreme accuracy. We set the 
working precision to several hundreds during the calculations to get some reliable numerical 



information for the model. If the effective potential is a parabola in Eq. (3.8) then the IR 
value of the dimensionless coupling is 

2 

w 



V = 1T 2 - (3-26) 
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If we take the parameters from Eq. ( |3,25 ), then this relation becomes the same. Let us 
note here, that the IR value ( |3.26| ) differs a factor of 2 of the one for which the pole occurs 
in Eq. ( |3.24 ). This is due to the circumstance that the relation (3.24) corresponds to the 
neglection of the higher-harmonics of the blocked periodic potential, whereas our numerical 
approach avoiding the Fourier-expansion takes automatically all higher harmonics with. In 
order to enlighten the deviation of our result from the Maxwell-cut induced parabolic 
effective potential (3.8) it is reasonable to plot 2 — yfi 2 as the function of the color N, see 
Fig.||. We succeeded to calculate it only up to N = 5. Nevertheless the figure clearly 
shows that y — > 2//3 2 for N — > oo. These results also give vanishing dimensionful coupling 




N 

Figure 3: The N dependence of the function 2 — y/3 2 . 



y in the IR limit. One can have the same conclusion as was obtained in the 3-color case, 
namely the low energy iV-color QCD2 is a massless interacting theory. 



4. Summary 

The phase structure of the bosonized version of the multi-color QCD2 is mapped out and 
it was shown, that the model possesses a single phase. 

After bosonization of the original fermionic model a periodic self-interacting scalar 
model is obtained with a mass matrix. The more involved last self-interaction term of the 
Hamiltonian ( |2.11| ) is shown to be treatable by Taylor-expansion giving further corrections 
to the mass matrix, so far the system is considered in or close to the ground state. The 
bosonized 2-color QCD2 coincides with the 2-flavor QED2 which implies the similar trivial 
phase structures of these models. The bosonized 3-color QCD2 and 3-flavor QED2 are, 
however, different models. The bosonized 3-flavor QED2 is known to be in the symmetry 
(periodicity) broken phase and represents a free massive theory. The scaling laws and the 
phase structure of the bosonized multi-color QCD2 has been determined by the functional 
RG technique, applying the Callan-Symanzik renormalization scheme. The periodic self- 
interaction has been found to be UV relevant, which tries to drive the flow into a pole, 
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signaling that the iV-color QCD2 is in the symmetry broken phase. The IR physics of the N- 
color QCD2 has been also determined after diagonalizing the mass matrix and integrating 
out the massive fields in the free-field approximation. It was found that the bosonized 
multi-color QCD2 represents effectively a SG-type model in the symmetry broken phase, 
characterized by a superuniversal dimensionless effective potential. In the IR limit the 
quantum fluctuations try to drive the system to a non-trivial saddle point what seems to 
be, however, never reached at finite energy scale. Nevertheless, the larger the number N 
of colors is the closer the dimensionless effective potential is driven to the parabolic shape 
J3.8|) . Making use of the bosonization relations, we have concluded that the iV-color QCD2 
is a massless interacting fermionic model. 
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